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We present a method for computing optical absorption spectra by means of a Bethe-Salpeter 
equation approach, which is based on a conserving linear response calculation for electron-hole co- 
herences in the presence of an external electromagnetic field. This procedure allows, in principle, for 
the determination of the electron-hole correlation function self-consistently with the corresponding 
single-particle Green function. We analyze the general approach for a "one-shot" calculation of 
the photoabsorption cross section of finite systems, and discuss the importance of scattering and 
dephasing contributions in this approach. We apply the method to the closed-shell clusters Na4, 
Na^ and Najj^, treating one active electron per Na atom. 
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I. INTRODUCTION 

Optical absorption spectra are a widely used and ver- 
satile tool for the characterization of the electronic prop- 
erties of systems of interacting electrons, including their 
elementary excitations in condensed-matter systems as 
well as in molecules and clusters. The information that 
can be obtained from absorption and transmission spec- 
tra complements the information accessible from photo- 
electron spectroscopy techniques, but the interpretation 
is more difficult because transitions in the interacting sys- 
tem are probed, and the correlated dynamics of electrons 
and holes in the many-electron system has to be taken 
into account in order to compare with experiments and 
to identify spectral characteristics. 

In a quasiparticle picture, the photoabsorption is de- 
termined by the probability of an electronic transition 
from an occupied to an empty (or virtual) quasiparticle 
state, and can be directly determined from the electron- 
hole correlation function. The determination of the ele- 
tron-hole correlation function is numerically demanding 
and usually effected by time-dependent density func- 
tional theory (TD-DFT) or by solving the Bethe-Salpeter 
(BS) integral equation derived in the framework of many- 
body perturbation theory^. 

In the latter approach, the BS equation is often used at 
the level of a random-phase approximation (RPA), and 
a popular choice for the calculation of the the excited 
states, which are needed as input for the BS equation, 
is the GW approximation for the self-energy, i.e., the 
effective interaction seen by one particle in the presence 
of the other carriers. (Here, G is the single-particle Green 
function and W the screened interaction.) 

In the case of state-of-the-art BS equation 
approaches^"™, which may be applied to extended 



and finite systems, the DFT single-particle states 
are used as input to GW calculations to obtain the 
quasiparticle correction to the Kohn-Sham eigenvalues. 
Then the equation for the electron-hole correlation is 
transformed into an eigenvalue problem, from which the 
real parts of the electron-hole transition energies are 
obtained. The numerical calculations needed for these 
Bethe-Salpeter approaches are simplified if one assumes 
an instantaneous, i.e., statically screened, electron- hole 
interaction for the direct part of the integral kernel 
of the BS equation. This approach works well with 
quite a few extended systemes, such as semiconductors 
and insulators. The straightforward inclusion of the 
full dynamic screening in the BS equation complicates 
its numerical solution tremendously, but it is possible 
to perform an expansion of the dynamically screened 
potential in plasmonic modcs^. Dynamic screening is 
also responsible for excitonic effects even in metals, such 
as Cu and Ag^. General attacks on the problem of 
dynamic screening in the BS equation go back to the 
Shindo approximation^. More recently this problem has 
been reexamined in the framework of nonequilibrium 
Green function theory, which allows one to derive a 
Dyson equation for the two-particle propagator—. 

In addition to standard BS equations there also exist 
explicitly time-dependent (nonequilibrium) Green func- 
tion based methods, which reduce to the BS equation at 
the GW-RPA level as a special case, but have been solved 
for simpler systems using more complicated two-particle 
correlations than those contained in the GW-RPA^ii^. 
It has also been demonstrated that Green function ap- 
proaches going beyond the GW-RPA accurately deter- 
mine the electronic energy renormalizations and broad- 
ening for collective excitations in extended system a^^'^^ . 

In finite systems, electronic states and electron-hole 
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transitions contributing to the absorption spectra are 
discrete, and are explicitly treated as such in most 
Green function based theories as well as configuration- 
interaction (CI) and TD-DFT calculations. However, 
the experimental signatures due to the discrete states 
and transition energies invariably appear with a fi- 
nite linewidth in photoemission and absorption measure- 
ments, respectively, due to the finite detector energy- 
resolution. In particular, the bunching of levels and tran- 
sitions in certain energy regions, which occurs with in- 
creasing system size, cannot be resolved experimentally 
and may lead to structures that are much broader than 
the experimental resolution, so that the underlying con- 
tributions from individual states or transitions cannot 
be distinguished. This is often accounted for by adding a 
phenomenological broadening, which is usually taken to 
be the same for all transitions. 

The approach of the present paper to the ab-initio cal- 
culation of absorption spectra in finite systems uses from 
the outset a finite broadening for single-particle states, 
i.e., finite lifetimes, as it is rigorously only the case in ex- 
tended systems. In addition to the background broaden- 
ing, however, there result dynamical interaction-induced 
features in the absorption spectra, namely linewidths 
and energy shifts, which are closely connected!^. We 
aim here at the consistent inclusion of both transition- 
energy renormalizations and finite linewidths due to ele- 
tron-hole correlations for finite systems. We achieve this 
by deriving an integral equation from a quantum ki- 
netic description of the dynamic electron-hole coherences 
driven by an external, coherent optical field, including 
Hartree-Fock (HF) and scattering contributions. Our 
approach is equivalent to an equation of the BS type 
for the complex electron-hole correlation function using 
the HF single-particle energies together with a frequency- 
dependent integral kernel, which contains correlation ef- 
fects beyond the static HF contributionsii. The finite 
linewidths in our absorption spectra exceed the "back- 
ground" single-particle broadening and thus describe in- 
teraction induced effects including the formation of new 
peaks. The broadening in our approach should be dis- 
tinguished from the artificial broadening in nonequilib- 
rium dynamicsiS. Since we base our derivation on the 
dynamics of electron-hole coherencesii in the interact- 
ing system, one can include self-energy effects and thus 
make the calculation self-consistent. An important prop- 
erty of the present approach is that it is straightforward 
to incorporate the requirement that the (two-particle) 
electron-hole correlation function be consistent with the 
electronic single-particle properties, so that a conserving 
approximation (in the sense of Baym and Kadanoff) re- 
sults. This approach fulfils important sum rules by con- 
struction^*. For instance, the f-sum rule for the density- 
density correlation, which can be read as the Thomas- 
Rciche-Kuhn sum rule for the absorption cross section in 
finite systems, is directly checked in our numerical calcu- 
lations. Last, but not least, we point out that including 
correlation effects in a consistent way by examining the 



equations for electron-hole coherences makes the interde- 
pendence between one- and two-particle properties clear. 
This approach therefore avoids possible inconsistencies 
of a two-step procedure of calculating first the single- 
particle properties, and then solving a BS equation with 
an effective two-particle interaction, as it is done in state- 
of-the-art BS equation approaches. 

Our approach to the ab-initio calculation of absorption 
spectra in finite systems is based on a recent development 
of a "consistent" GW theory for finite systems"'^^, and it 
has been applied to extract structural information from 
cluster absorption spectr In the present paper, we 
are mainly concerned with a detailed discussion of the 
general method for computing absorption spectra. We 
restrict ourselves to the derivation of a "one-shot" calcu- 
lational method using Hartree-Fock energies and a con- 
stant "background" broadening. This is appealing for its 
simplicity, and also because it preserves the distinction 
between HF and correlation contributions throughout the 
calculations. The present paper expands on the short 
discussion given in Ref J^, and a few illustrative applica- 
tions to small metal clusters, such as Na4, Na;^ and Najj^, 
which exhibit rather complicated correlation mechanisms 
between the interacting electrons and which develop true 
collective excitations with increasing clusters size. 

II. THEORY 

In this Section wc provide the derivation of the equa- 
tion for the electron-hole correlation function, which de- 
termines the photoabsorption cross section. We start by 
summarizing the connection between these quantities. 
We then provide the definitions and important proper- 
ties of the electron-hole correlation functions in Sec. Ill Al 
using non-equilibrium Green function theory. Sec. Ill Bl 
contains the derivation of the dynamical equation for the 
electron-hole coherence in the nonequilibrium formalism, 
and in Sec. Ill Ci this result is translated back into a two- 
particle equation for the electron-hole correlation func- 
tion. 

The central quantity for the calculation of the opti- 
cal absorption is the retarded density-density correlation 
function 

x'-(f,t;r',t') = (p(r,t)p(r',t')) - (p(f, t)) (p(r-, t')) (1) 

where p{r,t) = ip\r,t)ip{r,t) is the particle density op- 
erator, expressed by creation and annihilation operators 
■0^ and V'^i respectively. We use the expansion 

X'{r,t;?,t')= ^ (nin3|x'(i,i')l"2n4) (2) 

ni . . .n4 

^Vn, {r)'fn2 {rj^n, {T^)fn^ (^) 

in a basis of molecular orbitals {(/3„(r)}, where n la- 
bels the orbital. The density-density correlation also de- 
scribes the linear response of the system to weak pertur- 
bation, which is described here by the coupling of a field. 
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with matrix element in the orbital basis t/„j„2(^)j to the 
system Hamiltonian i?sys via H{t) — i?sys + ^^ext(i) with 



(3) 



The matrix element occurring in Eq. ^ can be expressed 
as the functional derivative 



(n2ni|x''(*7^')l"-3"-4) 



(4) 



for the determination of Calculating the kinetic and 
spectral Green functions including their dependence on 
two times is numerically extremely demanding, so that 
it is only feasible for simple systems^ii^. Instead of such 
a two-time calculation we would like to determine the 
optical response of a system at T = K or in equilibrium 
in terms of a quasiparticle spectrum^ i.e., we would like 
introduce an approximation that replaces the two-time 
Green functions by 
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where PnmaC^) — (^)cn2 (^) ^-^d c„ (cj^) denotes the 
annihilation (creation) operator for an electron in the or- 
bital ifn- Since the external field couples only electron 
and hole states, and in linear response the electronic dis- 
tribution functions p„„ are unchanged, the only nonzero 
elements in Eq. are those with indices that pair occu- 
pied (electron) and unoccupied (hole) states, i.e., Eq. dU 
is the electron-hole correlation function. The photoab- 
sorption cross section thus takes the formic 



G<itl,t2) 



G<{t) = G<{t,t) 

G^h-t2) 

1 



(8) 
(9) 
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G'{cj)e''"^'^-*'^duj 



1 CJ ^ 



In other words, in such an approximation the time depen- 
dence ti — t2 essential for the quasiparticle spectrum, is 
separated from the time dependence t = {ti+t2)/2 essen- 
tial for the kinetics of the quasiparticles. This is achieved 
by the generalized Kadanoff-Baym ansatz, which reads 
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ni . . .714 



Im{n2n,\x'\u:)\n3n^)idn,n, ■ dn,n,), G^n^i^M) ^ ^f^G'n,n,{h - t2) Gf,^^^{t2) 



(10) 



where 



(5) 



(6) 



ihG"^.. 



,.(h)Gl^,Jh-t2) . 



are the electric dipole matrix elements in the molecular 
orbital representation. 



A. Density-density correlation function 

The problem in the determination of the absorption 
cross section is the calculation of x'' {'^)j i-^-i the Fourier 
transformation in t — t' of Eq. (j?]). To obtain directly an 
equation for x^i^)j we take the detour of deriving dynam- 
ical equations for the density response, which are capa- 
ble of describing more general non-equilibrium situations, 
and adapt these for the response to weak perturbations. 
This approach is already implicit in the original treat- 
ment of conserving approximations for the two-particle 
correlation functionM, and has recently been exploited 
numericallyiiii^. The functional derivative in the defi- 
nition of the response function in Eq. Q can, in prin- 
ciple, be carried out by solving a dynamical equation 
for the density (pnn'it)) in the presence of an exter- 
nal field Unn'{t')- The non-equilibrium dynamics can be 
calculated in any many-particle formalism; here we use 
nonequilibrium Green function a^^'^^ ', because they offer 
the choice to introduce renormalized quasiparticle prop- 
erties at the single-particle level. For arbitrary nonequi- 
librium situations, this formalism yields coupled dynam- 
ical equations for kinetic and spectral Green functions, 
i.e., G^{ti,t2) and G''(ti, ^2), which depend on two times. 



We then end up with a kinetic equation for G^{t) that 
includes the quasiparticle properties via its dependence 
on G'^(w), i.e., the Fourier transform of G" with respect 
to ti — t2- 

In our recent paper, Refj^^, we were interested in the 
determination of the electronic quasiparticle properties, 
i.e., renormalized electronic energies and electronic life- 
times, in the GW approximation, and used different ap- 
proximations for the density-density correlation function 
X- Here, we summarize the connection of the quasiparti- 
cle properties with x ^nd the screened Coulomb poten- 
tial. The quasiparticle properties are directly determined 
by the interacting Green function 



[Go-i(l,2)-I](l,2)] G(2,l')-<5(1-1'), 



(11) 



where a space-time integration, or equivalently, a sum- 
mation over the single-particle indices and a time in- 
tegration, for quantities with an overbar is understood. 
Eq. (jlip depends on the non-interacting Green function 
Go (here taken to be the HF Green function) and the self 
energy (here taken to be the GW self energy) 



£(1,2) =ihG{l,2)W{2,l) 



(12) 



(In the orbital representation, Eq. is given by 

Eq. (fT9)) below.) In terms of the density-density correla- 
tion function, which is equal to the reducible polarization 
function, the screened interation W is an auxiliary quan- 
tity given by 



M^(l,l')=«(l,l')+«(l,2)x(2,3)«(3,l') 



(13) 



even though we need only the time-diagonal 
{p{t)) = -zhG<{t,t) 



(7) 



Note that the quasiparticle properties depend on the two- 
particle correlation function x via Eqs. pT|) - P^ . while 
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the two-particle correlation function depends on the 
single-particle properties via the generalized Kadanoff- 
Baym ansatz, so that the single- and two-particle func- 
tions are connected and the full problem becomes in prin- 
ciple self-consistent. 

The following remark on the role of the orbital repre- 
sentation used above is in order. First of all, since the 
dipole matrix elements in the HF representation connect 
electron states, i.e., unoccupied states above the LUMO, 
with hole states, i.e. occupied states below the HOMO, 
the quantities {pnmiit)) in Eq. Q describe coherent 
electron-hole amplitudes in the presence of an external 
potential U. We will refer to these expectation values, 
which are nonzero only in non-equilibrium, as electron- 
hole "coherences." The definition of these coherences de- 
pend on the orbital representation. Here, we always use 
the representation in HF molecular orbitals. This will 
prove particularly important in the next section, because 
it allows us to conveniently separate the HF contribution 
in the equation for the electron-hole coherence. 



B. Kinetic equation for the density response 

Our general approach to the calculation of the Green 
function G is discussed in some detail in Ref.^^ Thus we 
concentrate here on the derivation of the density-density 
correlation function. To calculate the action of the weak 
external field U, we start from the equation of motion 
for G^(ti,t2) evaluated at equal times ti = t2 = t which 
reads 

*ft^G<„,(i)= (14) 

"3 

~ ^ninai^T^) [2n2n2'^"2n3 + Un3n2{t) + Sn3n2 (^)] | 
+ Snin2 {t)- 

Here, T is the non-interacting part of the Hamiltonian 
Hsys ^ T + V, and includes the core potential in our 
case. The scattering contribution describing correlations 
beyond the mean-field (HF) approximation is given by 

Sn.n, (i) = E / [ ^"in3 ^^n, {I *) (15) 

"3 

+G<„3(t,t)S>„^(f,t)-(^o^)]. 
Also, the instantaneous HF self-energy 

S^fn2(*l) ^-*^E<"l"3|^^|«2n4)G<„Jti,ti) (16) 

"3 "4 

+ ih^ (nin2|w|n3n4) G^^^^{ti,ti) 



has been separated from the dynamic correlation contri- 
bution via the identity 

j:tn.(tut2)=^n!n2itl)Sih~h) (17) 
± e (±(ti - t2)) [S> „^ {h,t2) - S< „^ (ti, i2)] . 

In Eq. (jl6p and below, v is the bare Coulomb matrix 
element, with the index structure defined as 

{nin2\v\n3ni) ^ y^VidVa <^ (riX^ (ra) (18) 
X v{\ri - f2\)(pnAn)VnAr2)■ 
Eq. ([T7|) is quite general and holds regardless of the level 
of approximation for the self-energy, which determines 
the functional form of E^. Due to the separation be- 
tween the HF and the dynamic correlation contributions 
to the self-energy, and since we are not attempting a self- 
consistent calculation, it is natural to use for the basis set 
functions {ip7i{r)} the HF molecular orbitals, as opposed 
to DFT-based approaches that lump the exchange (Fock) 
and correlation part together, and separate out only the 
direct (Hartree) contribution to the self-energy. 

Specifically, we use the GW approximation for the dy- 
namic correlation contribution 

ri3ri4 

(19) 

where matrix elements of the screened Coulomb interac- 
tion W have been introduced in analogy with Eq. (|18l) . 
Note that Eq. is just Eq. ([T^ written in the basis of 
the HF orbitals. 

For the two-time Green functions from the correlation 
term we employ the generalized Kadanoff-Baym quasi- 
particle ansatz (ITO|) . Note that we have not specified yet 
which G'^{t — t) we use to describe the spectral properties 
of the electrons. If one uses the results of a GW calcu- 
lation, the calculation becomes self consistent, because 
the screened interaction depends on x- We do not in- 
tend to pursue the self-consistent problem further in the 
present paper, so that we neglect the correlation correc- 
tions to the HF single-particle properties in the follow- 
ing. With this approximation the retarded and advanced 
Green functions at the HF level are given by 

»ftG^(HF)(^ _ ^-)e^p|_^g^^^(^ „ (20) 

For clusters, the corrections to the ground-state ener- 
gies are expected to be small so that a Green function 
constructed from HF eigenvalues constitutes a reason- 
able starting point. A non-zero quasiparticle broadening 
7 ensures the proper behavior of the HF retarded and 
advanced Green functions, and we use here and in the 
following the notation e„ = e^^ — 17. With this ap- 
proximation, Eq. (jl4p is a kinetic equation for the den- 
sity response, which has a parametric dependence of the 
electronic quasiparticle properties via Eqs. p^ - (P(I| and 
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which is closed in the sense that it does not depend on 
Green functions with off-diagonal time arguments. 

Before we can actually compute as the functional 
derivative of Eq. (fH)) we need a few more steps. The 
functional derivative essentially captures the linear de- 
pendence of the density on the external potential U. 
To first order in U , or, equivalently, for a weak exter- 
nal potential, only density fiuctuations, i.e., {pnn'{t)) = 
—ihG^j^,{t,t) with n / n' , are infiuenced by the field, 
whereas the distribution functions /„(i) = {Pnn{t)) are 
unchanged and therefore time-independent. Specializing 
to the case of such a weak U we can rewrite Eq. (|14p as 

(z^^ - + eJ^f)G< „,(t) -t- (/„, - /„j[^C/„,„,W 



{{nin3\v\n2n4) - {nin2\v\n3n4))G^^^Jt) 



(21) 



Here, we have introduced the HF eigenvalues e^^^ and dis- 
tributions /„ = /(ej^^) where fifujj) denotes the Fermi 
function, so that /„ is equal to 1 for occupied and for 
unoccupied states at T = OK. In Eq. (PT|) . S stands for 
the scattering term of Eq. ((T5|) . which includes now the 
generalized Kadanoff-Baym ansatz. With the assump- 
tion of a weak external field, we should also consistently 
linearize the scattering term with respect to the G^_^^, 
and replace factors of , t) by their equilibrium value, 
ihf{e^^). With these modifications, the linearized ki- 
netic Eq. (pij) for the density response assumes a form, 
for which we can compute the functional derivative. 



C. Bethe-Salpeter equation for x"^ 

The equation for x'^(t, t') is obtained by functional dif- 
ferentiation of Eq. (pij) with respect to U{t') and letting 
[/ — ^ afterwards. This is done by replacing everywhere 
the term <5G< „^ (t)/<5[/„3nJt') with -ih{nin2\x''{t - 
i') 171371,4), cf. Eq. (HI). Apart from the SG/6U, functional 
differentiation of Eq. (PT|) also yields contributions of the 
form dW/SU but these are neglected, in agreement with 
other GW based approaches- The result depends only 
on ti —t2, and we can Fourier transform in the time dif- 
ference to obtain the frequency-dependent x'X'^) i^i the 
form 



with the correlation kernel 

{nin2\A{uj)\n3n4) 

r /ni(l ~ fns)fne + (1 ~ /rii)/ii5(l ~ fne) 



HF , 



^HF 



^HF 



cHF 



ni ' ' ^{ini + 7n4 ) 

X {ninz\v\nzn(i){niin2\v\nr^ni) 

fn^i^ ~ fn5)fne + (1 ~ fn2)fns{^ ~ fne) 



hjjJ 



pHF 



fHF 



^ni7i 



X {nin5\v\n3ne){nen2\v\n5n4) 

fnji^ ~ fne)fn5 + (1 ^ fnyjfnei^ ^ fn^) 



^ ^-^ huj — e 



HF , 



cHF 



cHF 



X {n4n5\v\nrnQ) {nen2\v\n5nj) 

fnri^ ~ fne) f 11;, + (1 ^ fn-r)fne{^ ^ fn^) 



.HF 



e^+e^^f +*(7n.+7n4) 



X {nin5\v\nrne) (nenslvlmm) 



(23) 



The functional form of Eq. ((22|) for allows one to 
draw conclusions how many-particle effects will show up 
in the calculated absorption spectra. Note first that the 
direct and exchange Coulomb contributions renormal- 
ize the transition energies together with the real part 
of the correlation (or scattering) contribution A. The 
imaginary part of A describes an interaction-induced de- 
phasing of electron- hole coherences driven by U. Thus 
the imaginary part of A determines the line broaden- 
ing of the absorption resonances. Note that the complex 
quantity A contains contributions to resonance shifts as 
well as the broadening; both effects are therefore related. 
The renormalization/broadening contributions A are de- 
termined by a coupling of each electron-hole transition 
to all other transitions due to the Coulomb interaction. 
Although the single-particle broadening 7 (inverse life- 
time) of each electron or hole state, or more generally, 
the single-particle spectral functions, influence the broad- 
ening of the transitions, Eq. ([25)1 shows that this influ- 
ence is not described by a simple convolution of the elec- 
tron and hole lifetimes. Rather, the line broadening for 
all absorption resonances arises due to an interaction- 
induced coupling of individual electron-hole transitions. 
The "collective" origin of the damping should be distin- 
guished from the question whether a given resonance is 
due to collective electron-hole excitations, such as ex- 
citons or plasmons, or whether it mainly stems from a 
single electron-hole transition. 



\n1n2\x'' {LL>)\n3n4:) 



(22) 



~^ (/"l fn2 ) 



+ 



(^{nin5\v\n2n6) - {nin2\v\n5ne)^ {n^nelx''' {i^)\n3ni) 

n^ne 

+ ^ (nin2|A(w)|n5n6)("-5?i6|x' ('^)l»^3»^4) =0 



III. RESULTS 

In this Section, we apply Eqs. ([5]) and ([^ together 
with Eq. ([25]) to compute the photoabsorption cross sec- 
tions for the clusters Na4 , Na^ and Na and we compare 
our calculations with experimental results. The pho- 
toabsorption cross section of Na4, Na^ and Najj^ has 
been calculated by different approaches, including: the 



jellium model approachl^"^ , highly correlated quantum 
chemistry methods such as multireference (double) Cl2^, 
a combination of GW and BS equation approaches'^, 
time-dependent density-functional theory (TD-DFT) in 
the local density approximation (LDA)^^ and beyond 
LDA with sophisticated exchange-correlation function- 
als2S. Another approach, which has been applied to cal- 
culate the photoabsorption cross section of small sodium 
cations (including Na,]"), is the TD-LDA (including gen- 
eralized gradient corrections) used along phase-space 
trajectories obtained within finite-temperature Born- 
Oppenheimer local spin density molecular dynamics'^. 

We focus on the influence of correlation effects due to 
the Coulomb interaction between the electrons from the 
outermost electronic shell of the Na ions at T =0K, and 
employ the frozen-phonon approximation, so that the ge- 
ometrical structure of the clusters is fixed at the ener- 
getically optimized configuration. Both core-shell elec- 
trons and electrons from the outermost shell influence 
the many-particle states of the interacting electron sys- 
tem and thus the optical spectra. However, because the 
core electron states are energetically well separated from 
the outermost shell, these states are only weakly mixed 
into the electronic states contributing to transitions in 
the energy window of interest, which is a few eVs, and 
we expect that the core electrons only introduce an ad- 
ditional fine structure in the optical spectrum in that 
window. For the clusters investigated here, we therefore 
treat one active electon per Na atom and use effective 
core potentials to describe the contribution of the core 
electrons. If one is interested in accounting for processes 
involving electrons from inner shells a two-particle-hole 
Tamm-Dancoff approximation'^*' or an algebraic diagram- 
matic construction^ may be used to describe ionization 
and double ionization processes in the inner valence re- 
gion. 

We treat here only closed-shell systems although, in 
principle, the theory is applicable to any finite system as 
the index n refers to a molecular orbital that can be from 
an open shell system as well. Since we are interested in 
the absorption lineshapes and interaction induced shifts 
due to electronic correlation effects, we view the struc- 
ture optimization of the cluster geomtry, and the nec- 
essary quantum-chemistry only as a tool. A relatively 
simple approach that fits well with our method is to ob- 
tain the HP eigenfunctions and the optimization of the 
underlying cluster geometries from a restricted HF cal- 
culatioii^i performed with the GAUSSIAN 03 quantum 
chemistry package^S-. In our calculations, the 2>s^ valence 
electron of each Na atom is explicitly treated using the 
lanl2dz basis set: a double zeta basis set of Gaussian- 
like atomic orbitals in the (3s3p/2s2p) contraction and 
relativistic Los Alamos^^ effective core potential (ECP) 
for the core electrons. The structure optimization of the 
clusters starts from geometry configurations taken from 
Refs^^>29^. 

For the computation of the correlation contribution A 
from Eq. (1^5)1 we use 7=0.25 eV for the quasiparticle en- 
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FIG. 1. Measured (a) and computed (b) absorption cross sec- 
tion for Na4. The experimental result is adapted from Ref.''— . 
The theoretical results are obtained with (solid line) and with- 
out (dotted line) correlation contributions in the equation 
for "x^ . The spectrum obtained using only the real part (plus a 
constant imaginary part) of the correlation term is also shown 
(dashed line). The inset shows the cluster geometry with dis- 
tances in A. 

ergies e ~ e — ij. The value for 7 used here is slightly 
modified from our earlier GW results, where we obtained 
for energetically low occupied states 0.085 eV for Najj" 
and 0.29eV for Najj^^^. However, this choice simplifies 
the numerical calculation because it makes it easier to 
carry out the sums over complicated energy denomina- 
tors, and makes A{uj) a smoother function of uj. We 
stress that introducing a fixed quasiparticle broadening 
for the calculation of A.{uj) does not mean that we are fix- 
ing the broadening of the resonances in to this value. 
Rather, the resonance broadening, which is responsible 
for the finite width of the peaks in the absorption spec- 
trum, is due to the imaginary part and the complicated 
index dependence of the full correlation kernel A(w). 

Figure [T] shows the computed absorption cross sec- 
tion for Na4 together with the experimental result from 
Refj2^. The main measured peaks which occur around 1.8 
and 2.5 eV are well resolved by the theoretical spectrum 
when x'^ is calculated by Eqs. ((22|) and ((23)) . To empha- 
size the importance of the correlation contribution, we 
compare this with the cross section calculated using a 
constant, phenomenological broadening for each transi- 
tion. Technically, this is achieved by using Eq. ([22)) with 
a constant, diagonal broadening, i.e., 

(nin2|A(w)|n5n6) = iS (24) 

where we have assumed F — 0.03 eV. Using the phe- 
nomenological broadening, the absorption cross section 
does not compare well with the experimental data as re- 
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FIG. 2. Measured (a) and calculated absorption cross sec- 
tions for NaJ. The experimental results are adapted from 
Ref.^. The spectra (b)-(d) are calculated for six different 
cluster geometries. The insets show the cluster geometries 
with distances in A. 



gards the position and the shape of the lines. Only the 
inclusion of the full correlation term yields good agree- 
ment between the calculated and the measured spectra. 
This result illustrates that the shift of the absorption res- 
onances due to the real part of correlation contributions 
can be on the order of several tenths of eVs and this 
shift is needed to obtain agreement with experiment. In 
addition to the shift, the imaginary part of A, which 
is responsible for the broadening, redistributes spectral 
weight between the different peaks and leads to an over- 
all very different shape of the absorption spectrum for 
the two cases. To disentangle the effect of the broaden- 
ing on the redistribution of the spectral weight, we have 
performed a calculation with only the real part of A, 
plus a diagonal broadening. This yields a splitting of 
the peak at 1.8 eV into two new peaks. Also, the peak 
at 2.5 eV is blue-shifted and now contains most of the 
spectral weight. This is in agreement with results for 
extended system, where it is known that the broadening 
and level shifts in this type of complex BS equation are 
closely connected and cannot be separated by taking into 
account the real or imaginary part only. 



When comparing experimental and theoretical results 
one should keep in mind that the geometry configura- 
tion, which we obtain by structural optimization, is im- 
portant for the calculation, but not well known for the 
clusters studied in the experiment. Moreover, finite tem- 
perature and inhomogeneous broadening effects in ex- 
periments wash out the intrinsic spectral features. To 
account for the finite temperature effects and for the 
unknown cluster geometries which might occur in ex- 
periment, we have considered six different configurations 
for Na^J" with low ground state energy: three structures 
with symmetry and three with C2v symmetry, see 
Fig. [2] The total ground state energies of the clusters dif- 
fer by less than 0.3 eV. The comparison between the mea- 
sured (from Ref. — ) and the calculated spectra is shown 
in Fig. [21 The experimental spectra exhibit a splitting 
of the main absorption line into a larger peak centered 
around 2.68 eV and a smaller peak around 2.98 eV. For 
all the clusters considered here, the main absorption lines 
occur in the same energy interval as the measured lines, 
but the shape of the peaks does not compare as well 
with experiment as in the case of Na4. Both the C4v and 
the C2v clusters from panel (b) of Fig. [5] yield a single 
large peak at 2.74 eV. The C4v cluster from panel (c) 
reproduces best the experiment with respect to the po- 
sition and the shape of the absorption lines. Also, the 
C2v cluster from panel (d) of Fig. [5] is capable of resolv- 
ing the smaller experimental peak at 2.98 eV, but the 
larger experimental peak at 2.68 eV is not present in the 
theoretical spectrum. In an experiment, clusters with 
different configurations may contribute to the shape of 
the absorption spectrum, but the weight of the contribu- 
tions due to different possible configurations cannot be 
determined from experiment nor from a T = K the- 
ory. Given the potential problems for experiment-theory 
comparisons, we find a good overall agreement between 
the measured and the calculated spectra of Na^ . 

Figure IIIII shows the comparison between the exper- 
imental and the theoretical cross section of Na^^^. The 
measured spectrum displays a large absorption line cen- 
tered around 2.74 eV—. Two different cluster geometries 
were used in the calculation: one corresponding to a pro- 
late (i.e., elongated) cluster and one structure with Cgv 
symmetry. The qualitative agreement between theory 
and experiment is very good with respect to the posi- 
tion and the shape of the peak for the prolate Na^j^ . For 
the cluster with Cgv symmetry the agreement is not that 
good, although the main peaks are situated in the exper- 
imental 2.3-3.2 eV energy interval. 

The different quality of the agreement between the- 
ory and experiment for the Na4 and Na^j^ clusters is also 
an indication that cluster size plays an important role: 
For small clusters, the geometrical configuration is rel- 
atively well defined, and the photoabsorption cross sec- 
tion is characterized by a few discrete resonances. In the 
other extreme, when the system contains a large num- 
ber of atoms there are many more overlapping individual 
electron-hole transitions, so that the resonances are not 
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FIG. 3. Measured (a) and computed (b) absorption cross sec- 
tion for NaJ^. The experimental result is adapted from Ref.''^. 
The theoretical spectra are calculated for the prolate cluster 
(solid line), and the structure with Cgv symmetry (dashed 
line). The insets show the cluster geometries with distances 
in A. 



well resolved. Moreover, there may be different geometri- 
cal configurations with similar total energies, which may 
contribute to the experimental spectrum. For Nag", we 
are neither in the limit of a small nor a larger cluster so 
that there is a competition between the two trends: the 
spectral density of the single electron-hole transitions is 
not yet large enough as to yield featureless broadened 
peaks, but the number of the possible configurations with 
different electronic properties is large enough so that so 
that the cluster geometry present in the experiment is 
difficult to predict. 

Finally we wish to point out that we have checked nu- 



merically that all the computed spectra shown fulfill the 
f-sum rule for finite systems to better than 99% 

— duJUJ idnin2 ' dnsm ) (25) 

ni . . .n4 

X Im(?^2^^l|x''(^)l»^3»^4) ^^e, 

m 

where m is the electron mass and Ne is the total num- 
ber of electrons in the system. Eq. ([25]) is the Thomas- 
Reiche-Kuhn sum rule for the absorption cross section 
and is a measure for the overall quality of the absorption 
spectrum. 

IV. CONCLUSIONS 

We have presented an approach to compute the ab- 
sorption spectra for finite systems based on a linear re- 
sponse theory for the dynamical electron-hole coherences 
in the presence of an external field. Using a quasi- 
particle ansatz in the quantum kinetic equation for the 
electron-hole coherence, which includes HF and scatter- 
ing contributions, allows us to derive a BS equation for 
the electron-hole correlation function with a complex, 
frequency-dependent kernel. The latter is determined 
by the scattering (or correlation) contributions to the 
electron- hole coherence dynamics. One thus obtains ab- 
sorption spectra that include level-shifts of transition en- 
ergies and a broadening of the absorption peaks due to 
electronic interactions. The correlation induced shifts 
and the broadening are closely related since they are the 
result of the complex, frequency dependent kernel, and 
the broadening mimics the bunching of discrete levels in 
finite systems. We apply the method to three important 
Na clusters: Na4 and the magic number clusters Na^ 
and NaJi^, and compare both the calculated resonance 
energies and the shape of the absorption lines with ex- 
perimental results. 

G.P. and W.H. acknowledge support through the Pri- 
ority Programme 1153 of the German Research Founda- 
tion. 
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